Nature Plants

Agrobacterium-mediated transformation is the dominant method for plant transgenesis, yet it frequently produces multi-copy, structurally complex T-DNA insertions associated with transgene silencing, unpredictable expression, and genome instability. Here, leveraging a high-throughput phenotypic reporter, we systematically dissect how T-DNA vector architecture, plasmid biology, and regulatory element choice shape transformation outcomes in Arabidopsis thaliana. We discover a pronounced trade-off between transformation efficiency and T-DNA copy number, uncovering the virulence enhancing overdrive sequence as a major determinant of this relationship. Guided by these insights, we engineered a new T-DNA vector that balances efficient transformation with predominantly single-copy integration. Additionally, we replaced viral elements, such as the widely used CaMV 35S promoter, with Arabidopsis-derived regulatory elements to minimise undesired enhancer effects, and developed a streamlined workflow for efficient T-DNA insertion mapping in the genome. Together, these advances form the T1 vector series, an Arabidopsis-optimised T-DNA vector system that enables clean, single-copy, and readily mappable transgene integration with predictable expression in the first generation after transformation.

Plants acquire iron using two canonical mechanisms: reduction-based uptake (Strategy I) and chelation-based uptake (Strategy II). Although Strategy II is characteristic of grasses, genes associated with chelated-metal transport occur broadly across plant lineages, obscuring how this pathway evolved. Here, we reconstruct the evolutionary history of Yellow Stripe-Like (YSL) transporters and Nicotianamine Synthase (NAS) across Archaeplastida and their closest non-plant homologs. Phylogenomic, distributional, and structural analyses reveal pronounced temporal uncoupling between these gene families. YSL most likely originated early in Viridiplantae and represents the deepest evolutionary module of the chelated-metal transport system. In contrast, NAS appeared much later through independent horizontal gene transfer events from fungi into euphyllophytes and specific moss lineages. This timing indicates that YSL-mediated transport initially functioned independently of nicotianamine, implying ancestral use of alternative siderophores. A pronounced expansion of YSL genes in Poaceae coincides with the emergence of canonical Strategy II, suggesting that this pathway arose by exaptation of pre-existing chelated-metal transport mechanisms.

Perennial crops follow different domestication trajectories from annuals, yet the molecular basis of slow-variable domestication--subtle tuning of conserved regulatory hubs--remains poorly characterized. We reconstructed the evolutionary history of the chloroplast kinase ABC1K7 across 9 seed plant species spanning [~]350 Myr, employing PAML codon models, IQ-TREE robust codon models, and protein-level phylogenetic inference, with AlphaFold2 structural modeling. ABC1K7 was under extreme purifying selection ({omega} = 0.073-0.104) across all seed plants. In coconut, a single Y[->]F substitution at residue 652--located >30 [A] from the catalytic core in a predicted intrinsically disordered region--represents the only non-synonymous change differentiating coconut from 7 of 8 angiosperm orthologs, and exhibits perfect co-segregation with domestication traits across a 17-year breeding panel (n = 327). These findings provide population-level evidence consistent with the slow-variable domestication model, identifying ABC1K orthologs as targets for perennial crop improvement.

Whole-genome duplications (WGDs) reshape plant genomes by generating redundancy, after which lineage-specific architectures emerge through fractionation, gene loss and rearrangement. How specialized metabolic pathways remain functionally integrated after such large-scale restructuring remains poorly understood. This problem is especially relevant for biosynthetic gene clusters (BGCs), which physically organize specialized-metabolism genes yet can be disrupted by post-duplication rearrangement. Here, we present the first chromosome-level genomes for Loganiaceae, including near telomere-to-telomere assemblies of Strychnos ignatii and S. pubescens, together with a draft genome of the extinct species S. ridleyi. Following a lineage-specific WGD, the two extant Strychnos species evolved contrasting genome-evolutionary trajectories and metabolite profiles: S. ignatii shows expansion of monoterpenoid- and monoterpene indole alkaloid (MIA)-associated gene families and strychnine-type MIA dominance, whereas S. pubescens exhibits elevated transposable element activity associated with DNA-binding with one finger (DOF)-linked regulatory rewiring and broader sesquiterpenoid- and triterpenoid-rich chemistry. Crucially, both species retain active strychnine biosynthesis despite fragmentation of a deeply conserved alkaloid BGC in MIA-producing Gentianales, revealing how pathway function can persist after disruption of ancestral BGC architecture. Comparative metabolomic and transcriptomic pathway analyses indicate norfluorocurarine oxidase (NO) as a major divergence point associated with strychnine accumulation. Promoter analyses, yeast one-hybrid assays, and electrophoretic mobility shift assays support a model in which S. ignatii retains the canonical jasmonate-responsive MYB, MYC2/bHLH, and AP2/ERF cis-regulatory module at NO, whereas the orthologous S. pubescens promoter shows reduced capacity to recruit these activators and instead exhibits a DOF-associated architecture. Together, our results show that WGD can decouple physical cluster architecture from pathway function, allowing specialized metabolic pathways to remain active while divergent chemical phenotypes evolve through lineage-specific combinations of coding-space expansion and transposable-element-associated cis-regulatory rewiring.

To overcome the enzymatic limitations of Rubisco, algae operate CO2-concentrating mechanisms (CCMs) that deliver concentrated CO2 to Rubisco tightly packaged in a specialized microcompartment called a pyrenoid. Pyrenoids are globally important biomolecular condensates, but their convergent evolution means that their molecular composition and emergent architecture cannot be inferred across clades. Here we characterize the pyrenoid of the Trebouxiophyceae alga, Chlorella sorokiniana. Using cryo-electron tomography, we provide an architectural overview of the pyrenoid and visualize pyrenoid-specific protein complexes. Quantitative proteomics and Rubisco co-immunoprecipitation followed by mass spectrometry demonstrate that inorganic carbon delivery machinery is conserved across green algae but the pyrenoid structural components are not. In vitro reconstitution supports the role of two previously undescribed proteins, one in assembly of pyrenoid traversing thylakoids (putative matrix thylakoid tether; PMTT) and another in starch tethering to the Rubisco matrix (putative matrix starch tether; PMST). In Nicotiana benthamiana, PMTT localized to the thylakoid stromal lamellae and PMST to chloroplast starch granules. Our findings provide insights into the molecular logic of pyrenoid assembly; how proteins mediate condensate-membrane and condensate-starch interactions; and expands the pyrenoid plant engineering toolkit, setting the stage for engineering a Chlorella pyrenoid into plants.

Cell polarity is an ancient organizing principle across kingdoms. As in animal epithelial cells, plant cells asymmetrically distribute proteins to establish functionally distinct membrane domains. In roots, radial polarity distinguishes inner and outer cell surfaces and supports directional nutrient transport, yet its molecular basis remains poorly understood. Here, we show that the leucine-rich repeat receptor-like kinases CaMRLK and IRK occupy complementary lateral plasma membrane domains in Arabidopsis thaliana roots. Polarity-guided proximity labeling identified previously uncharacterized proteins associated with inner- and outer-lateral domains. Clade VII LRR-RLKs, protein S-acyltransferases, SICK, IRKI1, and a distinct group of NPH3/RPT2-LIKEs assemble into the Lateral Protein Cluster (LPC) through multivalent interactions. LPC components are conserved across land plants, and disruption of NRL function impairs morphogenesis in Arabidopsis and Marchantia polymorpha. Together, these findings establish the LPC as an evolutionarily conserved molecular machinery linking radial cell polarity to plant morphogenesis.

Reported CRISPR/SpCas9 editing efficiencies in crops span 0-100%, but no quantitative synthesis has separated taxonomic family from delivery method, ploidy, clustering or publication bias. This meta-analysis estimated pooled per-T0-line editing efficiency across Cucurbitaceae, Brassicaceae, Solanaceae and Poaceae, and tested whether family is an independent moderator after adjustment for delivery and ploidy. A PRISMA 2020 systematic review identified peer-reviewed studies using SpCas9 with extractable per-line T0 edit counts; data were extracted independently by two reviewers, with inter-rater agreement reported. Logit proportions were synthesised with a binomial-normal generalised linear mixed model, and the family-as-moderator hypothesis was tested by a small-sample CR2 cluster-robust F-test on a three-level model with study-level clustering. Publication bias was assessed by Eggers regression and trim-and-fill. Twenty-two studies contributed 172 per-line effect sizes (Cucurbitaceae k=14, Brassicaceae k=20, Solanaceae k=68, Poaceae k=70). Pooled editing efficiency was 61.8% (95% CI 54.5-68.6%) with I{superscript 2}=93.4% ({tau}{superscript 2}=3.21) and a 95% prediction interval of approximately 5-98%. Per-family estimates ranged from 47.8% (Poaceae) to 73.8% (Brassicaceae); the univariate Q test was significant (p=0.0016), but family did not survive cluster-robust adjustment (F=0.73, p=0.63). Intraclass correlation placed 64.4% of variance at the study level, and Solanaceae remained dominated by a single study (58/68 rows). Funnel asymmetry was severe (Egger p<0.0001), and trim-and-fill reduced the bias-adjusted estimate to 45.2% (95% CI 39.0-51.5%). Apparent crop-family differences dissolve once within-study clustering and methodological covariates are accounted for; the bias-adjusted pooled estimate is closer to 45% than to 62%, and reported editing efficiencies reflect study-level factors more than taxonomic family. Key MessageApparent between-family differences in CRISPR/SpCas9 editing efficiency across four crop families reflect within-study clustering and publication bias, not intrinsic biology; family is not an independent moderator after cluster-robust adjustment.

Highly diversified specialized metabolism enables plant communication with pollinators, herbivores, and protectors1-4. Its chemodiversity, under which evenness, richness and variation is summarized5-7, includes many compounds without known function1-4 and presents an evolutionary conundrum about how and what is selected for8,9. Due to its complexity, it is frequently unknown how it is encoded in genomes. To produce population level chemodiversity, the traits need to allow for highly chemodiverse and highly specific individuals in the same population. Here we use metabolomics, transcriptomics, and genomics combined with field analyses and functional assays of monoterpene synthases in the Asteraceae Tanacetum vulgare (tansy) and identify forces which produce high population level chemodiversity: selection for product specificity in enzymes, loss-of-expression alleles, absence variation, and specialized metabolism islands drive individuals towards low chemodiversity while unlinked enzyme loci, expression variation alleles, presence variation, and de novo enzyme evolution enable high individual chemodiversity. Since the molecular data suggests selection for mechanisms that increase chemodiversity itself at the population level, the screening hypothesis which posited plants produce a reservoir of diverse chemicals prior to selection8 should be replaced by a chemodiversity selection hypothesis. The results demonstrate that, in addition to plant protection via individual chemicals with known targeting mechanisms for predators, being different from your neighbors even if you are closely related is likely an important element in plant protection.

Plants have evolved unconventional shoot systems with a "fuzzy" character that blurs organ boundaries, but the developmental programs driving this morphological diversity remain a mystery. One-leaf plants of the genus Monophyllaea are a striking example of this diversity: they develop a single indeterminate cotyledon instead of a conventional shoot system to form a unique structural unit termed a phyllomorph. Here, we integrate tissue-specific and single-nucleus transcriptomics with spatial gene expression analyses to characterize the meristems underlying this developmental system. The Monophyllaea-specific meristem of the phyllomorph retains a conserved shoot apical meristem-like transcriptional core but has a leaf lamina program superimposed. This chimeric transcriptional state integrates normally distinct developmental programs alongside altered phytohormone regulation. The one-leaf morphology emerges through the rewiring of conserved meristem regulatory networks rather than the passive suppression of meristem activity. The findings highlight the remarkable plasticity of developmental modules in plant morphological evolution.

Auxin is a key phytohormone that regulates all aspects of plant growth, development, and environmental responses, making the precise analysis of its distribution and signaling essential for understanding plant adaptation and physiological processes. However, despite the agricultural importance of oilseed rape (Brassica napus), the lack of robust, species-specific molecular tools limits detailed studies of hormone signaling in this crop. Here, we developed and characterized reporter systems for the sensitive visualization and quantification of auxin distribution and signaling in B. napus. The DR5cc auxin signaling reporter and a novel synthetic auxin-responsive reporter, BIP3, assembled from promoter fragments of three oilseed rape IAA genes, were generated to drive GUS expression. In hairy roots, both reporters showed auxin-responsive expression in the root apical meristem that became broader after auxin treatment. In transgenic seedlings, flowers at anthesis, and 12-day-old embryos, DR5cc exhibited a more defined expression pattern than BIP3. To monitor real-time auxin dynamics under abiotic stress, DR5cc fluorescent reporters were employed in hairy roots. Mannitol and NaCl treatments induced a time-dependent increase in fluorescence, peaking at 6-12 h before returning to basal levels after 24 h. Furthermore, dual-reporter assays enabled simultaneous monitoring of auxin and cytokinin signaling, revealing distinct hormone-specific spatial responses in hairy roots. Finally, we established a quantitative DII (qDII) reporter system using degron domains from B. napus Aux/IAA proteins, providing a high-resolution quantitative readout of auxin depletion. Together, these reporter systems enable spatial, temporal, and quantitative analyses of auxin dynamics during development and stress adaptation in oilseed rape.

Vascular plants rapidly coordinate root and shoot responses to water stress. Abscisic acid (ABA) mediates these adaptations; however, it remains unclear which cells produce ABA, whether ABA synthesis shifts during stress, and whether ABA movement is required for its adaptive functions. Here, we map ABA biosynthesis at cellular resolution in Arabidopsis and report that water-stress adaptive responses in roots and shoots require movement of ABA and its precursor AB-aldehyde from vascular tissues to target cells. We suggest that ABA accumulation arises from two parallel routes: (i) ABA synthesized in the vasculature via ABA2 and AAO3, then moving to guard cells, and (ii) phloem-derived AB-aldehyde being converted to ABA in the epidermis or bundle sheath by AAO1 and AAO2. Finally, we predict that tightly packed cells beneath leaf veins facilitate efficient ABA delivery to guard cells, an anatomical arrangement that has enabled angiosperms to evolve the use of ABA to rapidly close stomata.

Plant biology now publishes thousands of experimental research articles each year, but their core research logic, namely what questions are being asked, with what methods, and what is being found, remains locked inside free text and invisible to systematic analysis. Here we present a structured, 20-year atlas of The Plant Cell in which every paper is converted into a typed, directed Research Process Graph (RPG) of Question (Q), Method (M) and Finding (F) nodes connected by Q[-&gt;]M and M[-&gt;]F edges. A benchmarked large language model pipeline applied to 2,633 Plant Cell research articles published 2005-2026 recovered >110,000 Q/M/F nodes and >126,000 directed Q[-&gt;]M[-&gt;]F chains with>98% precision. A second LLM pass generalises each node into a paper-independent canonical form and assigns it to one of 10 top-level (L1) and [~]90 sub-level (L2) categories for each node type, producing the first comprehensive map of plant-biology research logic at the resolution of individual research questions. The atlas reveals that Plant Cell papers fall into seven canonical paper recipes with characteristic Q[-&gt;]M[-&gt;]F sub-structures, that peripheral experimental techniques have largely turned over while a stable methodological core persisted, and that the strongest correlate of per-PI citation impact is methodological breadth, not productivity or topical breadth. We release the atlas as a public, browsable database with five complementary interfaces: paper views, an LLM-powered research assistant, expert profiles, a taxonomy browser, and a method explorer. The database, available at https://rpg.connectome.tools/, turns the literature into a queryable community resource.

A fundamental question in developmental biology is how highly complex, yet reproducible multicellular body plans form from a single cell. The multicellular haploid body of the land plant Marchantia polymorpha develops from a single isolated cell - the spore - that divides asymmetrically. This produces a small terminally differentiating basal cell, and a large proliferative apical cell that gives rise to the multicellular sporeling body. The genetic basis for the establishment of asymmetry, first within the spore, and later within the sporeling body has remained virtually unknown. Here, we show that the plant specific RHO-type GTPase - RHO OF PLANTS (ROP) - polarises to the spore basal pole to ensure spore division is highly and reproducibly asymmetric. We demonstrate this highly asymmetric division is necessary to specify the terminal differentiation of the basal cell. Furthermore, we show that ROP-mediated polarised outgrowth is required to establish asymmetry within the sporeling body derived from the apical cell. Our discovery highlights how ROP function confers developmental robustness and contributes to the establishment of asymmetry during plant development from a single isolated cell.

Decoding how plants integrate multiple hormone signals to coordinate growth requires tools capable of resolving pathway interactions at cellular resolution in living tissue. Here we present ACE (Auxin-Cytokinin-Ethylene) and ACE2, proof-of-concept single-locus reporters to simultaneously capture activity of multiple hormones. Deploying ACE alongside well-established reporters, exogenous hormone treatments, and reverse-genetic perturbations of hormone biosynthesis, signaling, and transport in three-day-old etiolated Arabidopsis seedlings, we dissect the spatiotemporal hierarchy governing primary root elongation and root apical meristem (RAM) size. We demonstrate that both ethylene- and cytokinin-triggered root growth inhibition involve a boost of TRYPTOPHAN AMINOTRANSFERASE OF ARABIDOPSIS1 (TAA1)-mediated auxin biosynthesis and AUXIN RESISTANT1 (AUX1)-dependent auxin redistribution. Two spatially distinct auxin responses underlie the respective root growth effects: ethylene expands TAA1-dependent auxin biosynthesis from the root vasculature into the epidermis and promotes AUX1-mediated auxin import into the transition and elongation zones to inhibit cell elongation, while cytokinin confines ethylene-dependent TAA1-boosted activity to the vasculature and drives auxin accumulation in lateral root cap cells to reduce RAM size. Together, these data establish a reciprocal regulatory loop between these hormones, positioning ethylene as a convergence node in auxin-cytokinin crosstalk, and cytokinin as a modulator of the ethylene-auxin interaction. Critically, the changes in cross-activated reporter patterns described for different genetic backgrounds, alongside quantitative assessment of hormone-specific inhibition of the mutants growth, were consistent with the multi-hormone network established over two decades of research, and added cell-type-resolved spatial detail and a proposed hierarchy for the etiolated seedling root. Finally, a second-generation reporter, ACE2, overcomes key technical limitations of ACE, expanding the platforms capacity toward a higher-order multi-hormone monitoring system. These resources expand the Arabidopsis genetic toolkit and provide a generalizable framework instrumental for dissecting multi-hormone signaling hierarchies at the cellular level.

Plants detect neighbours through a reduced red-to-far-red ratio (R:FR), triggering elongation growth that reduces crop yield. Although Gibberellin (GA) is required for the neighbour-proximity (NP) elongation response, bioactive GA levels do not increase sufficiently to account for elongation magnitude, suggesting GA sensitivity as an additional regulated variable. Here, we show that GID1C, one of three Arabidopsis GA receptors, is the primary GA receptor involved in NP-induced elongation. GID1C protein accumulates selectively in hypocotyls and root tips under low R:FR without an increase in bioactive GA. The gid1c mutant shows a reduced elongation response that exogenous GA treatment cannot rescue. Transcriptome profiling reveals that GID1C controls 86% of the NP-responsive transcriptome, including genes for cell growth, division, and transcriptional regulation. Hub analysis identifies ICE1 as a GID1C-repressed transcriptional brake. ICE1 transcript is suppressed under low R:FR in a GID1C-dependent manner, and a phosphorylation-resistant ICE1 allele blocks NP-induced elongation. Together, these findings establish GA perception as an additional regulatory layer in NP, with subfunctionalisation among GID1 paralogs shaping the response to neighbouring plants.

The Arabidopsis PWWP-DOMAIN INTERACTOR OF POLYCOMBS1 (PWO1) and Telomere Repeat-Binding Proteins 1-3 (TRB1-3, TRBs) associate with distinct and shared protein complexes involved in epigenetic regulation, yet their cooperative roles in chromatin control and plant development remain largely unexplored. Here, we show that the interaction between PWO1 and TRBs is evolutionarily conserved. Both PWO1 and TRBs associate with plant telomeres, interact at these regions, and are co-enriched at subsets of interspersed telo-box motifs across regulatory regions genome-wide. TRBs facilitate PWO1 binding at shared genomic regions, including telo-box motifs. PWO1 and TRBs share a substantial number of genomic targets and preferentially bind chromatin regions associated with transcriptionally active states, whereas TRBs alone associate with repressive marks at thousands of loci. Genetic analyses show that the pwo1 trb1 trb3 triple mutant displays severe developmental defects, including main stem arrest and early maturation associated with aberrant lignin deposition in interfascicular tissues. In the triple mutant, key enzymes in the lignin biosynthesis pathway are upregulated, indicating that PWO1, TRB1, and TRB3 cooperatively regulate secondary cell wall formation. Together, our findings provide new insights into how PWO1 and TRBs cooperate to regulate chromatin states and orchestrate plant development, highlighting their central role in controlling gene expression programs. Significance statementThis study shows that PWO1 and TRB proteins co-occupy telomeres, including interspersed telo-box motifs, to regulate chromatin organization and plant development, particularly ectopic lignin deposition. Our findings reveal how these nuclear protein factors coordinate epigenetic states in Arabidopsis thaliana, providing a framework for understanding the control of developmental programs. O_FIG O_LINKSMALLFIG WIDTH=198 HEIGHT=200 SRC="FIGDIR/small/740627v1_ufig1.gif" ALT="Figure 1"> View larger version (38K): org.highwire.dtl.DTLVardef@1370171org.highwire.dtl.DTLVardef@3faea5org.highwire.dtl.DTLVardef@e30921org.highwire.dtl.DTLVardef@16c948e_HPS_FORMAT_FIGEXP M_FIG O_FLOATNOGraphical abstract:C_FLOATNO Evolutionarily conserved PWO-TRB interactions and their shared roles in chromatin regulation and plant development. Created with BioRender.com. C_FIG

The phytohormone auxin is a central coordinator of plant growth and development. Besides its canonical effect on gene transcription1,2, auxin triggers an ultra-rapid calcium ion influx that initiates the root gravitropic response3. The nature of the so-called rapid auxin pathway connecting the AFB1 auxin receptor3,4 and plasma membrane calcium channels remained unknown. Here, we show that auxin induces the direct interaction of the AFB1 receptor with the CNGC14 calcium channel. As the AFB1 receptor is independent of the ubiquitin ligase complex5, the auxin-induced interaction translates into relocalization of the receptor to the plasma membrane. We identify the interaction interface and provide evidence that the docking of the receptor to the channel complex activates Ca2+ influx and triggers growth inhibition. These findings position a calcium channel as an unprecedented component of the AFB1 auxin receptor complex. The ligand-dependent localization shift of a TIR1/AFB family receptor represents a novel paradigm in signal transduction and opens the possibility of unforeseen branches of auxin signaling pathways.

Nitrogen (N) is a major determinant of plant growth and productivity. Because soil N availability and internal N demand fluctuate, plants have evolved sophisticated mechanisms to coordinate N acquisition and utilization at the whole-plant level. However, how this coordination is achieved remains poorly understood. Here, we show that N-inducible LATERAL ORGAN BOUNDARIES DOMAIN transcription factors LBD37, LBD38, and LBD39 (LBDs) function as repressors of local N uptake and assimilation and systemic N-demand signaling in Arabidopsis. Triple mutants lacking these three LBDs displayed enhanced nitrate influx and increased accumulation of nitrate, amino acids, and total N. Transcriptome analysis identified an array of N-starvation- and nitrate-inducible genes derepressed in shoots and roots, including C-TERMINALLY ENCODED PEPTIDE (CEP) and CEP DOWNSTREAM (CEPD) genes, as well as genes involved in N uptake and assimilation. Grafting and genetic analyses revealed that LBDs gate the systemic N-demand signaling relay by repressing CEP and CEPD expression organ-autonomously. We also found that LBDs locally repress genes involved in N uptake and assimilation through a distinct regulatory mechanism. We propose that LBDs are key transcriptional repressors in a regulatory framework for optimizing N acquisition and utilization under fluctuating N conditions at the whole-plant level.

The domestication of maize from teosinte involved dramatic remodeling of the ear, yet the cellular and genetic bases of this transformation remain unclear. Here, we generate a single-nucleus and spatial transcriptome atlas of developing maize and teosinte ears. Comparative analysis reveals divergence in cob-associated cell types, with enhanced cytokinin signaling and reduced growth-inhibitory signals collectively driving cob thickening and enlargement in maize. We further demonstrate that domestication expanded the spatial expression domain of key transcription factors in maize meristem cells, enhancing the potential for increasing kernel number. Additionally, we verified a major domestication gene, ZmSPD1, in which two nonsynonymous SNPs differentiate maize from teosinte and alter jasmonic acid (JA) levels in the ear, thereby suppressing spikelet abortion to effectively double kernel production. These findings provide a cell-resolved mechanistic framework for how cob architecture and kernel number were shaped during maize domestication, offering new insights into the formation of key agronomic traits.

Hydathodes are specialized leaf structures present across vascular plants that allow guttation by connecting the xylem to the external environment through epithem tissue and permanently open water pores. However, the genetic mechanisms controlling their formation and physiological roles remain poorly understood. Here, we identify a genetic regulatory network that controls hydathode formation and links this process to leaf morphogenesis. This network converges on auxin signaling to coordinate the formation of the three hydathode cell types. We further show that epithem development requires sustained cell proliferation with limited endoreduplication. Analysis of hydathode mutants demonstrates that hydathode size and number are required to prevent reversible leaf flooding. Together, these findings establish a genetic framework for hydathode morphogenesis and uncover a central role for hydathodes in maintaining leaf water homeostasis under fluctuating environmental conditions.